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Abstract 

In  this  paper  we  construct  an  explicit  constitutive  model  that  is  capable  of  describing  the 
thermomecbanical  response  of  a  shape-memory  alloy.  The  model  consists  of  a  Helmholtz 
free-energy  function,  a  kinetic  relation  and  a  nucleation  aiterion.  The  &ee-energy  is  associated 
with  a  three-well  potential  energy  function;  the  kinetic  relation  is  based  on  thermal  activation 
theory;  nucleation  is  assumed  to  occur  at  a  critical  value  of  the  appropriate  energy  barrier.  The 
predictions  of  the  model  in  various  quasi-static  thermomechanical  loadings  are  examined  and 
compared  with  experimental  observations. 


1.  Introduction.  Thermoelasticity  theory  has  been  used  to  study  certain  general  issues 
pertaining  to  solids  that  undergo  reversible  stress-  and  temperature-induced  phase  transitions.  In 
this  theory,  the  potential  energy  function  that  characterizes  the  material,  as  a  function  of 
deformation  gradient,  has  multiple  energy-wells  for  certain  ranges  of  stress  and  temperature,  and 
each  energy-well  b  identified  with  a  phase  (or  variant  of  a  phase)  of  the  material.  For  example, 
the  studies  reported  in  [1, 8, 12, 19, 21, 22,  31, 34]  address  issues  pertaining  to  stable  equUibrium 
configurations  of  such  materials,  while  [2-S,  16, 18, 35]  examine  questions  related  to  evolution. 

A  complete  constitutive  theory  that  describes  the  behavior  of  such  materials  consists  of  a 
Helmholtz  free-energy  function  which  describes  the  response  of  each  individual  phase,  a 
nucleation  criterion  which  signals  the  conditions  under  which  the  transition  from  one  phase  to 
another  commences,  and  a  kinetic  law  which  characterizes  the  rate  at  which  this  transition 
progresses.  Explicit  examples  of  these  ingredients  have  been  constructed  b>  various  authors.  For 
example,  Ericksen  [13]  and  Silling  [34]  have  constructed  three-dimensional  Helmholtz 
free-energy  functions  that  model  certain  crystals.  Falk  [14]  has  studied  a  one-dimensional 
polynomial  free-energy  function  and  Jiang  [20]  has  used  a  similar  characterization  in  anti-plane 
shear.  Models  of  kinetic  relations  have  been  constructed,  for  example,  by  Otsuka  et  al.  [30]  by 
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relating  phase  boundary  motion  to  dislocation  motion  and  by  MUller  and  WUmansky  [25]  by 
using  certain  statistical  considerations. 

Recently.  Abeyarame  and  Knowles  [6]  presented  simple  models  for  each  of  these 
ingredients:  their  Helmholtz  £ree-energy  function  was  associated  with  a  piecewise  linear 
material,  while  their  kinetic  relation  was  based  on  thermal  activation  theory;  they  took  nucleation 
to  occur  at  a  critical  value  of  “dx  iving  force".  They  also  made  a  qualitiative  comparison  of  their 
theory  with  certain  experiments  on  some  shape  memory  alloys. 

The  present  study  generalizes  the  model  of  Abeyarame  and  Knowles  [6].  The  potential 
energy  function  in  [6]  has  at  most  two  energy  wells,  and  thus  it  only  accounts  for  two  phases. 
Howevar,  even  in  the  simplest  one-dimensional  setting,  one  often  encounters  three  material 
configurations  —  a  parent  phase  and  two  variants  of  the  product  phase.  The  principal 
generalization  of  the  present  paper  is  to  construct  a  three-well  energy  function  analogous  to  the 
two-well  mode!  of  [6].  This  allows  us  to  simulate  some  experiments  such  as  the  one  carried  out 
by  Ebrenstein  (see  Acbenbach  and  MUller  [7]),  in  which  a  bar  composed  initially  of  equal 
amounts  of  two  martensitic  variants  was  subjected  to  a  slowly  oscillating  stress,  and  at  the  same 
time  was  heated  and  then  cooled.  At  various  stages  during  the  ensuing  process,  the  bar  involved 
two  martensitic  variants,  two  martensitic  variants  and  austenite,  one  martensitic  variant  and 
austenite,  and  pure  austenite.  A  three-well  energy  is  essential  for  modeling  such  a  phenomenon. 

As  in  [6],  the  kinetic  law  utilized  in  the  present  study  is  based  on  thermal  activation 
theory.  The  nucleation  criterion  that  we  adopt  here  is  based  on  a  critical  value  of  the  appropriate 
“energy  barrier"  and  is  motivated,  in  part,  by  the  models  used  to  describe  nucleation  in  the 
materiab  science  literature;  see  e.g.  Christian  [10],  Fine  [15]. 

In  Section  2  we  outline  the  version  of  thermoelasticity  that  we  intend  to  use.  In  Section  3 
we  construct  the  three-well  eneigy  function  and  discuss  the  stability  of  the  various  phases.  Next, 
in  Section  4  we  calculate  the  energy  barriers  associated  with  this  energy  function  and  use  them  to 
establish  the  nucleation  criterion  and  kinetic  relation.  In  Section  5  we  carry  out  a  number  of 
simulations  and  compare  them  with  experimental  observations. 
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2.  Preliminaries.  In  this  section  we  review  some  relevant  concepts  &om  the  continuum  theory 
of  thermomechanical  processes  within  a  purely  one-dimensional  setting  that  corresponds  to 
uniaxial  stress  in  a  bar.  The  bulk  behavior  of  a  thermoelastic  material  may  be  characterized  by  its 
Helmholtz  free-energy  per  unit  mas ;  v(Y.0).  where  y  is  strain  and  0  is  temperature;  the  stress  o 
jid  specific  entropy  T]  at  a  particle  are  then  constitutively  related  to  y.  6  through 

o  =  P¥y(Y.  0).  n  =  -Vq  (y,  0),  (2.1) 

where  p  denotes  the  mass  density  in  the  reference  configuration.  The  potential  energy  per  unit 
reference  volume  GCf.  0f  o)  of  the  material  is  deHneo  by 

G(r,  0.  o)  =  pv(y.  0)  -  oy,  (2.2) 

and  its  val*  ..  at  an  extremum  of  G(*;  0,  o)  coincides  with  the  Gibbs  free  energy  per  unit 
reference  volume  g<y,  0)  =  pv(y,0)  -  pvy(y,  6)y. 

In  order  to  model  a  material  that  can  undergo  a  thermoelastic  phase  transition,  the 
function  G<*;  0>  o)  should  have  multiple  local  minima  (“energy-wells”)  when  the  temperature 
and  stress  take  on  suitable  values;  the  corresponding  Helmholtz  free-energy  potential  i|r(a,  0)  will 
be  non-convex,  and  the  stress-strain  curve  characterized  by  o  =  pvy(y,  0)  will  be  non-monotonic. 
In  this  theory,  each  local  minimum  of  the  potential  energy  function  G,  and  therefore  each  branch 
with  positive  slope  of  the  stress-strain  curve,  is  identified  with  a  different  (metastable)  phase  of 
the  material. 

Suppose  that  G(«;  0,  o)  has  at  least  two  local  minima  corresponding  to  a  given 
temperature  0  and  stress  o,  and  let  y  =  y  =  ^0,  o)  and  y  =  f  (0,  o)  denote  the  values  of  strain 
at  these  two  enerf^y-wells.  Then  the  strain  at  a  particle  that  is  subjected  to  this  0  and  o  could  be 
either  y  or  ^  depending  on  which  energy-well  (i.e.  phase)  the  particle  is  in.  Let  x=s(t)  denote  the 
current  location  in  the  reference  configuration  of  a  particular  particle  of  a  bar;  suppose  that  the 
particle  immediately  to  the  left  of  x=s(t)  has  a  strain  y  while  the  strain  on  its  right  is  then  x  = 
s(t)  denotes  the  location  of  a  phase  boundary.  i.e.  an  interface  that  separates  two  distinct  phases 
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of  the  material.  During  a  quasi-st?'ic  process,  the  rate  of  entropy  production  r(t)  in  a  segment  of 
the  bar  that  contains  the  interface  x=s(t)  but  no  other  phase  boundaries,  and  which  is  at  a  uniform 
temperature  6(t),  can  be  shown  to  be 

r(t)  =  f(t)s(t)/0(t).  (2.3) 

where 

f  =  G(f;0,o)-  G(y;0,o);  (2.4) 

f  is  known  as  the  "driving  force”  acting  on  the  phase  boundary.  The  second  law  of 
thermodynamics  thus  requires  that  the  following  entropy  inequality  hold: 

fsaO.  (2.5) 

If  the  driving  force  f  happens  to  vanish,  one  speaks  of  the  states  (f ,  0)  and  (y,  0)  as  being  in 
‘'phase  equilibrium”  and  of  the  quasi-static  process  as  taking  place  “reversibly”.  If  G(t;  0,  o)  > 

G(Tt  0.  o),  then  f  is  positive  and  so  according  to  (2.5)  one  has  s  ^  0;  thus  if  the  phase  boundary 
propagates,  it  moves  into  the  positive  side  and  thereby  transforms  particles  from  (f ,  0)  to  (y,  0). 
In  this  sense,  the  material  prefers  the  smaller  minimum  of  G.  This  is  also  true  in  the  reverse  case 
when  (j(f :  0,  o)  <  G(T,  0.  o).  One  therefore  speaks  of  the  phase  associated  with  the  lowest 
energy-well  as  being  the  stable  phase. 

By  using  the  first  law  of  thermodynamics  one  can  show  that  the  heat  getierated  when  a 
unit  mass  of  material  changes  phase  from  (f ,  0)  to  (y,  0)  is  i/p  +  X  where  X  =  0  (^  - 1  i)  is  the 
latent  heat;  if  the  phase  change  occurs  under  conditions  of  phase  equilibrium,  then  f=0,  and  the 
heat  generated  is  X. 

Let  X  =  s(t)  denote  the  Lagrangian  location  of  a  phase  boundary  at  time  t.  As  particles 
aoss  this  interface,  they  transform  from  one  phase  to  another  at  a  rate  that  is  determined  by  the 
underlying  “kinetics”.  The  kinetics  of  the  transformation  control  the  rate  of  mass  flux  ps  across 


the  phase  boundary.  If  one  assumes  that  this  flux  depends  only  on  the  states  (Y.  B)  and  (^,  6)  on 
either  side  of  the  interface,  then  the  propagation  of  the  phase  boundary  is  governed  by  a  relation 
of  the  form 

s  =  v(Y,  t.0),  (2.6) 

where  the  kinetic  response  function  v  is  determined  by  the  material.  Alternatively,  since  the 
constitutive  relations  a  =  B)  and  o  =  ptj>y(Y.  B)  can  be  inverted  (separately)  for  each 

phase,  one  can  express  ^  and  y  in  terms  of  a  and  B,  and  thus  re-write  the  kinetic  law  (2.6)  in  the 


s  =  v(o,  B), 


where  the  function  v  depends  on  the  two  particular  phases  involved  in  the  transformation  and  is 

i 

different  for  each  pair  of  phases.  Finally,  by  substituting  the  inverted  stress-strain-temperature 

i 

relations  into  (2.4)  one  ckn  express  the  driving  force  acting  on  an  interface  between  a  given  pair 

A  I 

of  phases  as  f  =:  f  (B,  o);  this  in  turn  can  be  inverted  at  each  fixed  B,  and  so  the  kinetic  law  can  be 
expressed  in  the  form 


s  =  V(f,B).  (2.8; 

The  basic  principles  of  the  continuum  theory  do  not  provide  any  further  information  regarding 
the  kinetic  response  functions  V;  in  particular,  explicit  examples  of  V  must  be  supplied  by 
suitable  constitutive  modeling. 

A  particle  changes  its  phase  either  when  it  crosses  a  moving  phase  boundary  or  by  the 
alternative  mechanism  of  nucleation.  Consider,  for  example,  a  quasi-static  process  in  a  bar  which 
involves  only  a  single  phase  for  some  initial  interval  of  time,  and  two-phase  states  at  subsequent 
times.  Hie  kinetic  law  describes  the  evolution  of  e:dsting  phase  boundaries  and  therefore  is  only 
operational  once  the  bar  is  in  a  two-phase  state.  The  transition  of  the  bar  from  a  single-phase 
configuration  to  a  two-phase  configuration  is  controlled  by  a  "nucleation  criterion".  If  a  particle 
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is  to  change  phase  from  a  state  (Y>  6)  to  a  state  (f .  6)  through  nucleation,  some  critical  condition 
of  the  form  ^6,  o)  =  0  should  be  reached  at  that  particle,  where  $  depends  on  the  pair  of  phases 
involved  in  the  transistion.  Models  for  $  are  often  based  on  the  notion  of  ‘‘energy  baniers”,  as 
will  be  described  in  Section  4. 

We  close  this  section  with  a  brief  discussion  of  some  qualitative  features  of  the 
one-dimensional  model  that  is  to  be  constructed  in  the  sections  to  follow.  Consider  momentarily 
the  rAree-dimensional  continuum  theory,  and  suppose  for  purposes  of  discussion  that  the  material 
at  hand  exists  in  a  cubic  phase  (austenite)  and  an  orthorhombic  phase  (martensite);  an  example  of 
this  is  the  class  of  Cu-Al-Ni  shape-memory  alloys  studied  extensively  by  Otsuka  and  co-workers, 
e.g.  [28-30],  and  more  recently  by  Chu  and  James  [11].  In  view  of  certain  invariance 
requirements,  the  associated  three-dimensional  potential  energy  function  G  must  have  seven 
energy-wells  corresponding  to  the  austenite  phase  and  the  six  “variants”  of  martensite.  During  a 
uniaxial  test  of  a  suitably  oriented  single  crystal  specimen,  and  for  suitable  values  of  temperature 
6,  the  material  is  found  to  remain  in  the  cubic  phase  for  sufficiently  small  values  of  stress,  in  the 
orthorhombic  phase  with  the  long  side  of  the  crystal  parallel  to  the  tensile  axis  for  sufficiently 
large  tensile  stresses,  and  in  the  orthorhombic  phase  with  its  long  side  normal  to  the  tensile  axb 
(i.e.  in  a  different  variant  of  the  orthorhombic  phase)  fox  sufficiently  large  compressive  stresses. 
In  the  one-dimensional  theory  we  model  this  by  allowing  G  to  have  three  energy-wells  for 
suitable  6  and  o,  the  ones  at  the  largest  and  smallest  values  of  strain  correspond  to  the  two 
variants  of  martensite  just  desaibed,  while  the  one  at  the  intermediate  value  of  strain  is  identified 
with  austenite.  Since  the  variants  of  martensite  are  crystallographically  identical  to  each  other 
when  o  =  0,  the  energy-wells  corresponding  to  them  are  required  to  have  the  same  height  at  all 
temperatures  whenever  the  stress  vanishes.  Moreo  •,  all  three  energy-wells  should  have  the 
same  height  if  the  stress  vanishes  and  the  temperature  coincides  with  the  transformation 
temperature  6^.  At  higher  temperatures,  the  phase  with  greater  symmetry  (austenite)  is  usually 

preferred  over  the  low-symmetry  phase,  and  so  the  model  is  to  be  constructed  such  that  the 
austenite  energy-well  is  lower  than  the  martensite  wells  when  6  >  6j;  the  reverse  is  true  for  6  < 
The  crystallographic  similarity  between  the  variants  also  suggests  that  the  specific  entropy 
associated  with  them  should  be  identical,  and  therefore  that  the  latent  heat  associated  with  the 
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transformation  from  one  variant  to  another  should  be  zero;  this  too  is  a  feature  of  our  model. 

Finally  a  note  on  terminology:  for  simplicity  of  presentation  we  shall  sometimes  speak  of 
the  ‘‘three  phases”  rather  than  the  ‘‘one  phase  and  two  variants”;  similarly  we  shall  often  use 
the  term  ‘‘phase  boundary”  generically  to  refo'  to  both  an  interface  between  two  phases  and  to 
an  interface  between  two  variants  (which  ought  to  be  called  a  twin  boundary). 

3.  Helmholtz  Cree-energy  function.  In  this  section  we  construct  an  explicit  three-well 
Helmholtz  free-energy  function  \(r<Y.  6)  that  characterizes  the  response  of  a  multi-phase  material; 
the  three  energy-wells  are  viewed  as  corresponding  to  an  austenitic  phase  and  to  two  variants  of 
martensite.  Recall  fust  that  the  elastic  modulus  coefficient  of  thermal  expansion  a  and  specific 
heat  at  constant  strain  c  of  a  thermoelastic  material  are  related  to  \|r(Y,  6)  through, 

PVyy<Y.B)  =  H.  *V^<Y.0)/Vyy(Y.6)  =  a.  -ev00(Y.  6)  =  c.  (3.1) 

Thus  if  a  and  c  are  constant  on  some  domain  of  the  (Y,  6)-plane.  then  by  integrating  (3.1)  one 
finds  that 

PV<Y.  B)  =  <l«/2)  (  Y  -  g*  )  ^  -  1*0  Y  0  +  pc  (0  -  0  log(0/0*)J  +  pv*  (3.2) 

on  that  domain,  where  g*,  0«  and  are  constants. 

Consider  a  material  which  exists  in  a  high-temperature  phase  austenite  (A)  and  as  two 
variants  (M'*’  and  M')  of  a  low-temperature  phase  martensite.  Suppose  for  simplicity  that  the 
austenite  and  both  martensitic  variants  have  the  same  constant  elastic  modulus  (i>0,  the  same 
constant  coefficient  of  thermal  expansion  a>0  and  the  same  constant  specific  heat  oO.  (The 
model  that  follows  can  be  readily  generalized  to  describe  the  case  wherein  the  different  phases 
have  different  but  constant  material  properties.)  By  (3.2),  the  Helmholtz  free-energy  function 
t|r(Y,  0)  associated  with  this  material  must  have  the  form 


pV(y,0)  = 


<lV2)(Y-gi)^ 

<W2){Y-g2f 

<K®(T-83)^ 


jiayB  +  pc  ^8  -  0  log(0/0^)^  +PVi  o°  Pj* 
pa  Y  0  +  pc  ^0  -  0  log(G/02)^  +  PV2  0°  ^2  ’ 
paY0  +  ^0  ‘  0  log(0/02)^  +PV3  0°  P3' 


(3.3) 


where  p  is  the  mass  density  of  the  material  in  the  reference  cocfiguration,  and  0-,  g..  \|r.,  i=l,2,3. 
are  nine  additional  material  constants  whose  physical  significance  will  be  made  clear  in  what 
follows.  The  regions  P^,  P2  and  P^  of  the  <Y,  0)-plane  on  which  the  three  expressions  in  (3.3) 
hold  are  the  regions  on  which  the  respective  phases  A,  and  M'  exist;  they  are  assumed  to 
have  the  form  shown  in  Figure  1.  where  in  particular  the  boundaries  of  P^,  P2  and  P^  have  been 
taken  to  be  straight  lines.  The  temperature  levels  0^  and  Om  denote  two  critical  values  of 
temperature:  for  6>0m  the  material  only  exists  in  its  austenite  form,  whereas  for  0  <  0^  the 
material  only  exists  in  its  martensite  forms.  Throughout  this  paper  we  will  restrict  attention  to 
temperatures  less  than  V 

We  now  impose  a  number  of  restrictions  on  v  ^  oftlcr  that  it  properly  model  the 
stress-free  response  of  the  material  we  have  in  mind.  Since  the  potential  energy  function  G  and 
the  Helmholtz  energy  function  coincide  when  the  stress  vanishes,  any  characteristic  to  be 
assigned  tn  G  at  o  =  0  could  be  equivalently  imposed  on  \{r.  Let  0^  and  8m  be  the  two  critical 
temperatures  shown  in  Figure  1,  rnd  let  0^  denote  the  transformation  temperature,  0  <  <  Bj 

<  0j^.  We  assume  that  all  three  phases  A,  and  M’  may  exist  when  o  =  0, 0  =  0j.  Therefore 
the  function  \|K*.  0j)  must  have  three  local  minima,  with  the  minima  occuring  at  the  smallest, 
intermediate  and  largest  values  of  strain  identified  with  M',  A  and  respectively.  Since  0j  is 

the  transformation  temperature,  the  value  of  y  at  each  of  these  minima  should  be  the  same.  Next, 
for  o  s  0  and  0  close  to  0^,  we  require  y(*,  0)  to  continue  to  have  three  energy-wells.  Moreover, 

since  M***  and  M'  are  regarded  as  * ‘variants”  of  each  other,  the  two  martensite  snergy-wells  must 
have  the  same  height  for  all  temperature  in  this  range;  in  addition,  for  0  >  0^  the  martensite  wells 
must  be  higher  than  the  austenite  wrll,  while  for  0  <  0j  they  should  be  lower. 

On  enforcing  these  requirements  on  the  function  y  defmed  by  (3.3),  one  frnds  that 
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Vi-V2  =  Vi-V3"^^°* 


<3.4) 


log  Oyap  =  ”  <g2  ■  gp  +  '°8  (Bj/Bi)  =  ^  <S3  ■  gl>  + 

In  (3.5)  we  have  let  Xp  denote  the  common  value  of  \|f2  aod  ifTp-  tfr^;  one  can  readily  verify 
that  Xp  represents  the  latent  heat  of  the  austenite  •*  martensite  transitions  at  the  transformation 
temperature  &p  and  that  the  latent  heat  associated  with  the  M'*‘-  M'  transition  is  zero. 


The  stress-response  function  o<Y,  0)  =  PVy(Y.  0)  associated  with  (3.3)  is 


o(Y.  0)  = 


li<Y-gi) 

lt<Y-g2) 

^(Y-gj) 


-  pa0 

-  pa0 

-  |xa0 


on  Pp, 
on  P2. 
on  Pj. 


(3.6) 


Two  graphs  of  a(Y.  0)  versus  Y  are  shown  in  Figure  2:  Figure  2(a)  corresponds  to  a  fixed  value  of 
temperature  in  the  range  0j^<  0  <  6^  and  the  stress-strain  curve  shows  three  rising  branches 
corresponding  to  the  three  phases  A,  M'*’  and  M';  Figure  2(b)  is  associated  with  a  temperature  in 
the  range  0  <  0  <  0|q  atid  the  two  rising  branches  of  the  stress-strain  curve  correspond  to  the 
variants  M"*"  and  M'. 


Of  the  three  parameters  gp,  g2  and  g^,  one  is  fixed  by  the  choice  of  reference 
configuration,  while  the  other  two  are  determined  thrbugb  the  transformation  strains.  In 
particular,  if  we  choose  stress-free  austenite  at  the  transfiprmation  temperature  0^  to  be  the 
reference  state,  then  by  setting  c(0, 0j)  =  O  in  (3.6)p,  one  obtains 

gp=-a0p..  (3.7) 

' 

Next,  let  Y>p  (>0)  denote  the  transformation  strain  (see  Figure  2(a))  between  each  martensitic 
variant  and  austenite.  Then,  fiom  (3.6), 
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(3.8) 


Yt  =•- 82-81= 


(This  can  be  readily  generalized  to  the  case  where  the  transformation  strain  between  phases  A 
and  M'*'  is  say  Yp  and  that  between  A  .'md  M"  is  yj  »•  Yj  •)  Finally,  if  (3.4),  (3.5),  (3.7)  and  (3.9) 

are  substituted  back  into  (3.3),  one  finds  that  p\K(Y>  B)  contains  the  term 

+  (juJiya?"^  -f  pc6  log(6^/6j)  as  an  inessential  linear  function  of  temperature  which  may 

be  eliminated  by  taking 

pVj  =  -  (p/2)a2e|.  Bj  =  Op.  (3.9) 

In  summary,  the  material  at  hand  is  characterized  by  the  common  elastic  modulus  p, 
specific  heat  at  constant  strain  c  and  coefficient  of  thermal  expansion  a  of  the  phases;  the 
stress-free  transformation  temperature  6^;  the  mass  density  p  in  the  reference  state;  the  latent 

heat  Xj  at  the  temperature  0^;  and  the  transformation  strain  Yj.  The  Helmholtz  free-energy 
function  is  given  by  combining  (3.3)  with  (3.4),  (3.5),  (3.7)  ^  (3.9): 


PV(Y.0)  = 


(P/2)y^  -  paY(B-aj.)  +  pc  0 -  log(0/aj.))  on  Pj, 

(p/2)(  Y  -  Yj  )  ^  ’  pct(Y "  YxX0  -  0x)  +  pcti  (l  -  log(0/0.j.))  -  pX^(l  -  0/0.p)  on  P2, 
(|i/2)(  Y  +  Yj)  ^  -  pa(Y  +  YyXB  •  0^)  +  pc0  (l  -  log(0/aj.))  -  pX,j.(l  -  0/0,j.)  on  P3. 

(3.10) 


The  various  other  thermo-mechanical  characteristics  of  the  material  can  now  be  derived  from 
(3.10).  In  particular,  the  stress-response  function  o(Y,  0)  =  0)  is  given  by 


6(y.0)  = 


PY  -  pa(0-aj.)  on  P^, 

p(Y-Yx)  -  pa(0-0ij.)  on  P2, 
P(Y+Yx>  -  pa(0-aj.)  on  P3. 


(3.11) 


In  order  to  complete  the  description  of  the  Helmholtz  free-energy  function  we  need  to 
specify  the  regions  P^,  P2  and  P3  of  the  (Y.  0)-plane  shown  in  Figure  1,  i.e.  we  need  to  specify 
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the  boundary  curvec  y  =  Yj(0)  shown  in  the  figure.  To  this  end,  we  first  prescribe  the  stress-levels 

at  the  local  maxima  and  minima  of  the  stress-strain  curve.  As  indicated  in  Figure  2,  we  take,  for 
simplicity,  these  stress-levels  to  be  given  by  ±Oj^(0)  and  ±o^<0).  In  view  of  our  earlier 

assumption  that  the  boundaries  of  the  regions  P^,  P2.  P^  are  straight,  the  functions  and 

o„(0)  must  be  linear  in  0.  Moreover,  since  according  to  Figure  1  we  must  have  y9(0_)= 
V®m^’  **  =  0  and  Oj^(0j^j)  - 


Of/O)  =  liM(0  -  fOi  Bjq  <  e  <  Bj^, 

Oj^(B)  =  nm(0  -  Bj^)  +  mM(0j^  -  0^)  -  ny.j.  for  0  <  B  <  B^^, 


where  m  and  M  are  positive  material  constants.  The  equations  y  =  Y^(0)t  i=l. 2,3,4,  describing  the 
boundaries  of  P^^,  P2  and  P^  are  then  given  by  ±aj^(0)  =  o(  yj<0),  6),  i*3,2,  and  ±o^(ISi)  = 
o<ti<e).0),i=4.1. 

Thu?  far,  we  have  only  described  y  on  the  (“metastable”)  portion  ^1  +  ^2'*’  ^3 
(y>  B>plune.  It  is  not  necessary,  for  the  purposes  of  the  present  section,  to  specify  an  explicit 
form  for  y  on  the  remaining  (“unstable”)  portion  of  this  plane;  any  function  y  which  is  once 
continuously  differentiable,  is  such  that  is  negative  on  the  unshaded  portion  of  Figure  1 ,  and 

conforms  with  (3.3)  would  be  acceptable.  An  infinite  number  of  such  functions  exist,  provided 
only  that  the  material  parameters  satisfy  certain  inequalities;  this  is  discussed  in  the  appendix. 

Next,  it  is  useful  to  map  the  regions  P|,  i=l,2,3,  of  the  (y,  B)-plane  associated  with  the 
respective  phases  A,  and  M’  onto  the  (6,  o)-plane  using  o  =  o(y.  B).  The  result  of  this 
mapping  is  displayed  in  Figure  3  where  P-  is  the  image  of  P^.  Given  the  stress  o  and  temperature 
B  at  a  particle.  Figure  3  shows  all  of  the  phases  that  are  available  to  that  particle. 

Finally  we  address  the  issue  of  the  stability  of  the  phases.  The  potential  energy  function 
G(T>  B,  o)  =  p\|Ky»  B)  -  06  of  the  material  at  hand  can  be  calculated  using  (3.10).  At  each  (0,  o), 
G  has  one  or  more  local  minima.  Where  G  has  multiple  energy-wells,  one  can  use  the  explicit 


-11- 


formula  for  O  to  determine  the  particular  minimum  that  is  smallest.  This  determines  the  phase 
that  is  stable.  The  result  of  (his  calculation  is  displayed  in  Figure  4.  The  stress-level  o^(6) 
indicated  in  the  figure  u  given  by 


Oo®> 


5^1  /  e_ 

Y'j>  \ 


-  1 


(3.13) 


and  is  known  as  the  Maxwell  stress  for  the  A-M^  transition.  The  Maxwell  stress  for  the  A-M' 
and  M'^-M"  transitions  are  and  0  respectively.  The  two  states  A  and  M***  that  are  both 

associated  with  any  particular  point  on  the  boundary  o  »  both  have  the  same  value  of 

potential  energy  G  and  both  are  stable;  if  these  states  co-exist  and  are  separated  by  a  phase 
boundary',  the  driving  force  on  that  mterface  would  be  zero  and  the  phases  would  be  in  phase 
equilibrium. 

Suppose  momentarily  that  a  particle  always  chooses  the  phase  that  is  stable  from  among 
all  phases  available  to  it  Then  the  response  of  a  particle  as  the  stress  or  temperature  is  slowly 
varied  is  fully  determined  by  Figure  4.  For  example,  consider  a  fixed  temperature  6  greater  than 
the  transformation  temperature  6j.  As  the  stress  o  is  increased  monotonically  from  a 
sufficiently  negative  value,  the  particle  is  in  the  martensite  variant  M*  until  the  stress  reaches  the 
value  -Og(6):  it  then  transforms  to  austenite  and  remains  in  this  phase  as  the  stress  increases  over 
the  intermediate  range  <  o  <  o^<0);  at  o  s  o^<6)  the  particle  transforms  to  M***  and 

remains  there  for  o  >  o^<6). 

I  The  immediately  preceding  discussion  assumed  that  a  particle  is  always  in  the  stable 
piiase.  In  solids  however,  particles  can  often  remain  for  long  times  in  states  that  are  merely 
metastable  and  the  transition  from  a  metastable  phase  to  a  stable  phase  is  controlled  by  additional 
coluiderations,  viz.  nucleation  and  kinetics.  We  now  turn  to  these  issues. 
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4.  Nucleation  and  kinetics.  Given  the  stress  o  and  temperature  6  at  a  particle.  Figure  3  shows 
the  various  phases  that  that  particle  can  adopt,  while  Figure  4  indicates  which  among  them  is  the 
stable  phase.  If  a  particle  happens  to  be  in  a  phase  that  is  not  its  stable  phase,  the  questions  of 
whether,  and  how  fast,  it  will  transform  into  the  stable  phase  are  answered  by  a  nucleation 
criterion  and  a  kinetic  law.  In  this  section  we  will  describe  simple  models  for  nucleation  and 
kinetics  under  the  assumption  that  the  phase  transitions  are  "thermally  activated". 

A  particle  can  change  its  phase  in  one  of  two  ways:  a  particle  whose  phase  is  the  same  as 
that  of  its  neighboring  particles  could  change  its  phase  spontaneously  through  “nucleation”  if 
the  stress  and  temperature  at  that  particle  satisfy  an  appropriate  nucleation  criterion. 
Alternatively,  a  particle  in  one  phase  adjacent  to  a  particle  in  a  different  phase  and  separated 
from  it  by  a  phase  boundary,  will  change  phase  through  “growth”  if  the  phase  boundary 
propagates  towards  it,  the  motion  of  the  phase  boundary  being  controlled  by  a  kinetic  law. 

Energy  barriers:  Figure  5  shows  two  schematic  graphs  of  the  potential  energy  function 
G<Yi  6,  o)  plotted  versus  y  for  fixed  <0,  o).  Figure  5(a)  shows  three  local  minima  and 
corresponds  to  a  pair  (6,  o)  at  which  all  three  phases  co-exist  (i.e.  the  point  (0,  o)  lies  in  the 
region  common  to  P^,  ‘‘“d  P3  “  Figure  3);  Figure  5<b)  corresponds  to  a  pair  (0,  o)  at  which 

only  the  two  martensite  variants  co-exist  Figures  S(a)  and  S(b)  of  the  potential  energy  function 
correspond  to  the  respective  Figures  2(a)  and  2<b)  of  the  stress-strain  curve.  The  six  quantities 
by(0,  o)  indicated  in  Figure  5  are  the  energy  barriers  to  a  transformation  from  phase-  i  to 
phase-j,  where  it  is  convenient  to  use  the  subscripts  1, 2  and  3  to  refer  to  the  phases  A,  M'*’  and 
M'  respectively.  In  order  to  calculate  these  energy  barriers,  we  need  expressions  for  the  values  of 
G  at  the  local  maxima,  and  this  in  turn  requires  a  knowledge  of  the  Helmholtz  free  energy 
function  \{f  on  the  unshaded  regions  of  Figure  1.  While  we  can  suitably  continue  (3.10)  into  this 
“unstable  region”  in  many  different  ways,  our  present  purpose  is  merely  to  construct  a  simple 
continuum  model  that  describes  the  various  features  of  the  theory.  Consequently,  we  will  simply 
extend  the  parabolas  which  describe  G  near  its  local  minima  in  Figure  5  until  they  intersect,  and 
use  the  values  of  G  at  these  intersection  points  as  estimates  for  the  values  of  G  at  the  local 
maxima.  Using  (3.10)  to  calculate  G(Y5  0,  o)  =  p\jr(y,  0)  -  oy  and  then  carrying  out  this 
calculation  leads  to  the  following  expressions  for  the  six  energy-barriers: 
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b^3<0*  [  o  +  Ojj(0)  +  hYt/2  '?, 

b3j<e.o)  =  <2n)’^  [o  +  o^(e)-HYT/2  ]^. 

bj2<0. o)  =  <2^)’^  [o-o^(0)-M.Y-r/2 

b2i(0,o)  =  (2jJi)‘Mo.Ojj(0)  +  nYp/2  (4.1) 

b23(0.a)  =  <2n)'^  (o+  hYtI^ 

b32<9. o)  =  (2n)‘^  [o-  HYyl^ 

where  each  b|j(0.  o)  is  defined  for  values  of  <0.  o)  at  which  the  ith  and  phases  co-exist,  i.e. 
the  point  (0.  o)  lies  in  the  region  common  to  Pj  and  Pj  in  Figure  3;  0^(0)  is  the 
austenite-martensite  Maxwell  stress  introduced  previously  in  (3.13). 

Nndeation:  Considering  first  the  question  of  nucleation,  we  suppose  that  a  particle 
which  is  in  phase-i  will  transform  to  phase-j  by  nucleation  if  the  relevant  energy  barrier  bjj(0,  o) 
is  less  than  some  critical  number  n-jj  the  associated  nucleation  criterion  is  thus  given  by  setting 
bjj(0,  o)  =  n.j  in  (4.1).  In  view  of  the  symmetry  of  the  potential  energy  function  G  when  o  =  0, 

we  shall  assume  that  both  the  A-^M^  and  A-»M*  transitions  nucleate  at  the  same  value  of 

temperature  if  the  stress  vanishes;  a  similar  assumption  for  the  reverse  M'*’-»A  and  M*-^A 

transitions  will  also  be  made.  The  former  value  of  nucleation  temperature  is  denoted  by  M  (for 

s 

“martensite  start”)  while  the  latter  is  denoted  by  A  (for  “austenite  start”).  Finally,  we  will  also 

0 

assume  that  at  any  given  temperature  0,  the  nucleation  stress-level  for  the  transition  is 

the  negative  of  the  nucleation  stress-level  for  the  reverse  transition  at  that  same 

temperature.  One  can  readily  enforce  these  restrictions  on  the  n^j’s;  when  combined  with  (4.1), 
this  leads  to  the  following  nucleation  criteria  for  the  various  transitions: 

o  +  Oo(0)  forA-»M’, 

o  +  o^j(0)  ao^(A^)  forM'-»A, 
o-o^(0)  i-o^(M^)  forA-^M'*’, 

o-Oq(0)  s-Og(Ap  forM'^-*A,  (4.2) 
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o>2 


for 

o  s  -  2  for 

where  o  <0)  is  the  austenite-martensite  Maxwell  stress  given  by  (3.13)  and  the  constants  M  ,  A 

O  V  » 

and  2  are  characteristic  of  the  material;  necessarUy,  <  0j  <  and  2  >  0. 

Figure  6  Ulustrates  these  nucleation  criteria  on  the  (6,  o)-plane.  If  the  inequalities  in  (4.2) 
hold  with  equality,  the  resulting  equations  describe  a  set  of  straight-lines  in  the  (0,  o)-plane;  the 
nucleation  criterion  states  that,  as  indicated  in  the  figure,  crossing  one  of  these  lines  nucleates  an 
associated  transition.  Figure  6,  as  shown,  corresponds  to  a  material  for  which  the  critical 
nucleation  stress-level  given  by  (4.2)  for  the  M~-^A  transition  exceeds  the  corresponding 
stress-level  for  the  A-»M'''  transition  for  some  range  of  temperature,  i.e. 

2  >  (1/2)  (p>-j/y.j.aj.)  (A^  -Mp;  (4.3) 

this  need  not  necessarily  be  the  case. 

If  the  cunent  state  of  the  bar  were  to  involve  either  a  temperature  or  stress  gradient,  one 
could  determine  the  location  in  the  bar  at  which  nucleation  occurs.  In  this  paper  we  will 
consider  a  uniform  bar  that  is  always  subjected  to  uniform  stress  and  temperature  fields;  the 
location  of  the  nucleation  site  in  this  bar  is  therefore  rather  arbitrary.  If  the  bar  bad  been  rendered 
inhomogeneous  by  a  slight  uniform  taper  with  the  small  end  at  x=K)  and  the  large  end  at  x=L, 
then  the  transition  from  a  low-strain  phase  to  a  high-strain  phase  would  necessarily  commence  at 
x=0  and  the  reverse  transition  would  occur  at  x=L.  We  shall  arbitrarily  assume  that  this  is  the 
case  in  our  uniform  bar  as  well. 

Kinetics:  Let  x  =  s(t)  denote  the  location  at  time  t  of  a  phase  boundary  in  a  bar,  and  let 
(Y>  o)  and  (f,  0,  o)  denote  the  strain,  temperature  and  stress  at  the  two  particles  adjacent  to  the 
phase  boundary  on  its  left  and  right,  respectively.  Suppose  that  the  particle  on  the  left  is  in 
phase-i  while  the  particle  on  the  right  is  in  phase-j  (recall  that  1=1, 2, 3  corresponds  to  phases  A, 
M^,  M~  respectively).  The  driving  force  f  -  G(f ,  0,  o)  -  (}(Y,  0,  o)  on  an  i/j  phase  boundary  can 
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be  calculated  from  (4.1): 


f  =  Y.j,(  o  -  0^(0) )  foraM'VA  interface, 
f  =  -  Yj'(  o  -  0^(0) )  for  an  A/M'*’ interface, 

f  =  Yjf  o  +  0^(0)  )  for  an  A/M'  interface,  (4.4) 

f  =  -  Y<j'(  o  +  Oq(0)  )  foraM'/A  interface, 
fs2oYj  for  a  M'*^/M' interface, 
f=-2oYj  for  a M'/M"*^  interface. 

As  the  phase  boundary  propagates  through  the  bar,  the  particle  immediately  in  front  of  it 
"jumps"  from  one  local  minimum  of  G  to  another,  and  an  explicit  model  of  the  kinetic  relation 
may  be  constructed  by  viewing  this  jumping  process  on  an  atomic  scale.  In  order  to  jump  from 
one  minimum  of  G  to  the  other,  an  atom  must  acquire  an  energy  ai  least  as  great  as  that 
represented  by  the  relevant  energy  barrier:  for  an  atom  undergoing  a  phase-i-»^  phase-j  transition 
this  barrier  is  bs-(a,  0);  for  the  reverse  phase-j -»phase-i  transition  it  is  b-;(a,  0).  Under  suitable 

IJ  «  jj 

assumptions  about  the  statistics  of  this  process,  the  probability  that  the  energy  of  an  atom  is  at 
least  as  great  as  B  is  exp(-B/K0)  where  K  is  Boltzmann’s  constant.  The  average  rate  at  which 
atoms  jump  from  one  minimum  to  the  other  is  taken  to  be  proportional  to  the  probability  of 
exceeding  the  corresponding  energy  barrier;  we  assume  for  simplicity  that  the  proportionality 
factor  is  the  same  for  the  phase-i-»phase-j  and  phase-j-^phase-i  transitions.  The  velocity  s  of  the 
phase  boundary,  being  the  macroscopic  measure  of  the  net  rate  at  which  atoms  change  from 
phase-i  to  phase-j,  is  then  taken  to  be  the  difference  in  the  average  rates  associated  with  the  i->j 
and  the  j-»i  transitions: 


where  r  denotes  the  number  of  atoms  per  unit  reference  volume  and  R*j  is  a  positive 
proportionality  factor,  related  in  part  to  the  frequency  with  which  atoms  attempt  to  crossover  to 
the  new  phase. 
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Substituting  (4.1)  into  (4.5)  now  leads  to  an  explicit  representation  for  the  kinetic 
relations  of  the  various  transitions  in  the  form  s  =  V|j(o,  6).  By  using  (4.4),  they  can  be  expressed 
in  the  alternative  form  s  =  Vjj(f,  6) : 


2Rijexp  - 


^  liV/4  I 

J 


(4.6) 


where  g  =  Yj  for  both  M'*’/A  and  A/M"  interfaces,  and  g  =  2Y-p  for  an  interface.  These 

kinetic  relations  automatically  satisfy  the  condition  fVj^j(f,6)  >  0,  so  that  any  motion  consistent 
with  them  will  conform  with  the  entropy  inequality  (2.5).  Moreover,  at  each  fixed  6,  Vjj(f,6) 
increases  monotonically  with  f,  so  that  the  greater  the  driving  force,  the  faster  the  speed  of 
propagation.  If  the  driving  force  f  is  small,  so  that  quasi-static  processes  take  place  close  to  phase 
equilibrium,  one  can  approximate  (4.6)  to  obtain  the  following  kinetic  relation  which  is  linear  in 
f; 


s 


rKe 


exp 


I- 


l*g 

8rKe 


f. 


(4.7) 


5.  Thenno-mechanical  response  of  the  modeL  In  this  section  we  will  utilize  the  constitutive 
model  constructed  in  Sections  3  and  4  to  determine  the  uniaxial  response  of  a  bar  when  it  is 
subjected  to  various  thermo-mechanical  loadings.  We  describe  the  analysis  associated  with  one 
of  these  loadings  in  some  detail;  the  analysis  corresponding  to  the  remaining  cases  is 
conceptually  similar. 

Consider  the  isothermal  mechanical  loading  of  a  uniform  bar  at  a  temperature  0  >  A^. 

The  bar  is  initially  unstressed  and  is  composed  of  austenite.  As  the  stress  o(t)  is  mononotonically 
increased,  the  bar  remains  in  this  phase  for  some  time  0  <  t  <  t^.  By  (3.11),  the  elongation  5  of 

the  bar  during  this  stage  of  loading  (measured  from  the  reference  state)  is  given  by 


8(t)/L  =  o(t)/ji  +a(0-ftj.)  for  0<t<tj. 


(5.1) 


From  Figure  6  and  the  parag;raph  below  (4.3)  we  conclude  that  M'*'  martensite  is  nucleated  at  the 
lefit  end  of  the  bar  at  the  instant  t  =  t^,  where  t^  is  given  by 

o<tj)  =  o^<e)-o^(Mp.  (5.2) 

During  the  next  stage  t^  <  t  <  t2.  the  bar  is  composed  of  M'*'  martensite  on  the  interval  0  <  x  < 
s(t)  and  austenite  on  s(t)  <  x  <  L.  By  (3.11), 


8(0  =  s(t)  [  o(t)/n  +  Yj  +  a  (0  -  0.J,)  ]  +  t  L  -  s(t)]  [  o(t)/n  +  a  (0  -  ap)  ]  for  tp<  t  <  t2. 

(5.3) 

where  the  phase  boundary  location  s(t)  is  found  by  integrating  the  appropriate  kinetic  relation  in 
(4.6).  (4.4),  i.e. 


s(t)  2  R  exp 


f^(t)  +  li^Y^4  ] 

2\u^KQ  J 


s(tp)  «  0, 


(5.4) 


with  the  driving  force  f(t)  =  [a(t)  -  0^(0)]  Yj*  At  the  instant  t  =  t2.  the  phase  boundary  reaches 
the  right  end  of  the  bar,  s(t2)=L.  For  t  >  t2  the  bar  consists  entirely  of  M'*'martensite  and  its 
response,  according  to  (3.11),  is  given  by 

6(t)/L  a  o(t)/n  +  Yj  +  a  (0  -  0j)  fort>t2.  (5.5) 

A  similar  analysis  can  be  used  to  describe  the  response  corresponding  to  subsequent  unloading, 
as  well  as  to  loading  by  compressive  stress  in  which  case  the  M‘  variant  is  involved  instead  of 


In  order  to  study  the  quantitative,  as  well  as  qualitative,  suitability  of  our  model,  we 
chose  values  for  the  material  parameters  that  are  of  the  conect  order  of  magnitude  for  a 
Cu-Al-Ni  shape-memory  alloy: 
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(5.6) 


H  =  3  X  10^°  N/m^,  p  =  8000kg/m^ 

a  =  16  X  10'^/°K.  c  =  400JAg°K, 

y.p=  0.05,  Jtj.  =  5.7  X  10^  J/kg, 

Oj.  =  307®  K.  =  308®  K.  =  306®  K, 

2=^2.5xlo‘^N/m^ 

K  =  1.381  X  10*^^  J/®K, 

Rjj  =  0.0448  m/s, 

r  =  9.046  X 10^®  atoms/m^. 

The  values  of  the  transformation  and  nucleation  temperatures  6,|,,  A  and  M  were  taken  from 

IS  s 

Otsuka  et  al.  [281  and  correspond  to  an  alloy  whose  composition  is  Cu-14.0  Al-4.2  Ni  (wt%).  By 
comparison,  the  remaining  material  parameters  are  less  sensitive  to  the  alloy  composition  and  the  | 
values  chosen  for  them  do  not  correspond  to  an  alloy  of  this  precise  composition.  The  modulus  ^ 

I 

p,  the  transformation  strain  and  the  latent  heat  Xq.  at  the  transformation  temperature  were  | 
estimated  using  data  in  Otsuka  and  Wayman  [29].  The  values  of  the  mass  density  p,  coefGcient  j 
of  thermal  expansion  a,  and  rpecific  heat  c  were  estimated  using  information  in  [24].  The  value  | 
of  the  M  /M  nucleation  stress  Z  was  obtained  from  Sakamoto  et  al.  [32].  The  value  of  R^2 
estimated  by  using  our  kinetic  relation  (4.6)  in  conjunction  with  the  austenite/martensite  interface  ' 
velocities  measured  by  Grujicic  et  al.  [17]  and  reported  in  their  Figure  5;  the  remaining  R|j’s 
were  arbitrarily  assumed  to  have  this  same  value.  The  number  of  atoms  per  unit  reference 
volume,  r,  was  calculated  by  using  the  mass  density,  alloy  composition  and  the  atomic  masses  of 
Cu,AlandNL 

Tnere  are  four  other  material  parameters,  viz.  m,  M,  6^  and  V  that  are  involved  in  the 
description  of  our  model.  Even  though  they  do  not  affect  the  response  of  the  material,  the 
validity  of  the  model  requires  that  such  numbers  exist  in  a  manner  that  is  consistent  with  the 
various  constitutive  inequalities  given  in  the  appendix.  One  can  show  that  there  is  a  range  of 
acceptable  values  for  these  parameters. 

Results:  (i).  Figure  7  shows  two  force-elongation  curves  corresponding  to  isothermal 
mechanical  load  cycling  by  the  application  of  a  prescribed  stress.  Figure  7(a)  conesponds  to  a 
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temperature  below  M  ,  with  the  bar  transforming  between  the  M’  and  variants  without 

involving  austenite.  Figure  7<b)  corresponds  to  a  temperature  greater  than  A  ;  as  the  stress 

s 

increases  from  a  sufficiently  negative  value,  the  bar  transforms  from  M"  to  A  and  then  from  A  to 

+  *42 

M  .  The  loading  and  unloading  rate  underlying  both  of  these  figures  is  lo(t)l  =  5  x  10  N/m  sec. 

The  response  depicted  in  these  figures  is  similar  to  that  observed  by  Nakanishi  [27]  for  Au-Cd; 

see  Figures  10, 13  in  [27]. 

(ii).  Figure  8  shows  two  elongation-temperature  curves  which  result  from  cycling  the 
temperature  with  the  stress  held  fixed;  Figures  8(a)  and  (b)  correspond  respectively  to  tensile  and 
compressive  values  of  the  applied  stress.  The  bar  transforms  between  the  phases  A  and  M'*'  in  the 
former  case,  between  the  phases  A  and  M'  in  the  latter.  Observe  that  the  transformation  from  the 
high  temperature  phase  (austenite)  to  the  low  temperature  phase  involves  an  elongation  in  Figure 
8(a)  and  a  contraction  in  Figure  8(b).  The  heating  and  cooling  rate  underlying  both  of  these 

figures  is  IB(t)l  s  O.OOl^^sec.  The  response  depicted  in  these  figures  is  similar  to  the  response 
desaibed  by  Krishnan  et  al.  [23],  Figure  13  for  Cu-Zn  at  a  constant  tensile  stress  and  that 
observed  by  Burkart  and  Read  [10],  Figure  6  for  In-Tl  under  compressive  stress. 

If  the  initially  austenitic  bar  remains  stress-free  as  it  is  cooled  from  a  high  temperature, 
the  phases  and  M'  are  both  nucleated  simultaneously  at  6  =  M  ,  the  former  at  the  left  end  of 

the  bar,  the  latter  at  the  right  end;  they  then  grow  at  the  same  rate  according  to  their  kinetic 
relations,  and  once  the  transformation  is  complete,  the  bar  consists  of  an  equal  mixture  of 
and  NT.  Since  the  transformation  strains  involved  in  the  A->M^  and  A-*M'  transformations 
have  been  taken  to  be  and  -Yj  respectively,  the  length  of  the  bar  does  not  change  due  to 

transformation.  The  elongation-temperature  response  in  this  case  is  therefore  a  straight  line.  This 
is  a  trivial,  one-dimensional  manifestation  of  what  is  sometimes  referred  to  as  ‘*self 
accomodation”. 

(iii).  Figure  9  displays  the  result  of  a  calculation  which  attempts  to  model  the  “shape 

memory  effert”.  We  begin  with  a  martensitic  bar  which  is  composed  of  M"*"  for  0  <  x  <  L/2  and 

M'  for  1/2  <  X  <  L  and  whose  initial  temperature  is  less  than  M  .  The  bar  is  first  subjected  to  a 

s 
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program  of  isothermal  mechanical  loading  during  which  time  the  stress  is  first  increased  ^nd  then 

decreased  back  to  zero.  At  the  end  of  this  stage  of  loading,  the  bar  is  composed  of  only  (as 

can  be  deduced  from  Figure  6),  the  stress  has  returned  to  the  value  zero,  and  the  bar  has  suffered 

a  permanent  elongation.  During  the  next  stage,  the  bar  remains  unstressed  while  it  is  first  heated 

to  a  temperature  greater  than  A  (which,  according  to  Figure  6,  transforms  it  to  phase  A)  and  is 

s 

then  cooled  back  to  its  original  temperature  (which,  by  Figure  6,  leads  to  a  configuration 
involving  equal  amounts  of  and  M');  at  the  end  of  this  thermal  loading,  the  state  of  the  bar  is 
identical  to  its  original  state.  In  the  calculations  underlying  Figure  9  we  took  the  mechanical 

I* 

loading  and  unloading  rate  to  be  lb(t)l  =  8  x  10^  N/m^sec  and  the  heating  and  cooling  rate  to  be 

I0(t)l  =  0.08°K/sec.  Schematic  figures  similar  to  Figure  9  may  be  found,  for  example,  in 
Krishnan  et  al.  [23]. 

(iv).  Next  we  simulate  the  experiment  canied  out  by  Ehrenstein  as  described  by 

Achenbach  and  MUller[7].  Consider  a  martensitic  bar  which  is  initially  at  zero  stress  and  at  a 

temperature  <  M  ;  the  segment  0  <  x  <  L/2  of  the  bar  consists  of  M"  while  L/2  <  x  <  L  consists 
s 

of  M'*’.  Wc  consider  a  time  interval  0  s  t  ^  tp  and  apply  a  stress  o(t)  =  c ^(1  -  cos  2jit/80)  while 
simultaneously  varying  the  temperature  according  to  0(t)  =  0(0)  +  E2(  1  -  cos  2jit/tp).  We  take  0 

<  Ej  <  1/2  and  o(tp/2)  <  -  Oo(Ag)  which  ensures  that  the  hottest  temperature  0(tp/2) 

is  large  enough  to  nucleate  austenite.  The  resulting  elongation  history  is  shown  in  Figure  10(a) 
and  may  be  compared  with  Figure  2  in  Achenbach  and  MUller  [7];  the  loading  parameters 
underlying  our  figure  are  E^  =  5  x  10^  N/m^,  E2  =  4°K,  0(0)  =  304°K  and  tp  =  4000  secs. 

The  calculations  show  that  the  macroscopic  response  of  the  bar  plotted  in  Figure  l(Ka)  is 
associated  with  the  local  transitions  shown  in  Figure  10(b):  during  an  initial  stage,  roughly  0  <  t 

<  978  sec,  the  bar  consists  of  only  phases  M'  and  M***;  the  driving  force  on  the  M'/M***  interface 
is  negative  (by  (4.4))  and  this  causes  the  interface  to  move  leftward  in  accordance  with  the 
appropriate  kinetic  law  (4.6);  thus  during  this  stage  the  amount  of  M  ^  increases  at  the  expense  of 
M*.  For  the  heating  rate  used  in  our  calculation,  this  leftward  moving  interface  has  not  yet 
reached  the  end  x  =  0  of  the  bar  when  t  =  978  sec;  at  roughly  this  instant,  the  nucleation  criterion 
for  the  NT -»  A  transition  is  satisfied  and  phase  A  is  nucleated  at  the  left  end  of  the  bar.  As  t 
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continues  to  inaease,  the  newly  emerged  A/M'  interface  propagates  to  the  right  (since  its  driving 
force  is  positive)  in  accordance  with  its  kinetic  law,  while  the  M'/M'*'  interface  continues  to 
move  leftward;  the  amount  of  phase  M'  thus  continues  to  decrease  while  the  amounts  of  phases 
A  and  M***  increase.  A  short  while  later  (t  •  1026  sec)  the  nucleation  criterion  for  the  A 
transition  is  satisfied  and  phase  A  is  nucleated  at  the  right  end  of  the  bar.  There  are  three 
interfaces  in  the  bar  at  this  time,  viz.  a  rightward  moving  A/M'  interface  and  two  leftward 
moving  interfaces,  one  M'/M"*"  and  the  other  M^/A.  At  t  -  1136  sec  the  first  two  of  these 
interfaces  meet  so  that  during  the  next  stage,  the  bar  transforms  from  M^->  A  as  the  A/M^ 
interface  advances  towards  the  M'*'/A  interface.  Eventually  these  two  phase  boundaries  meet  at  t 
••  1314  sec  and  the  entire  bar  consists  of  phase  A.  During  the  next  stage  1314  <  t  <  3304,  the  bar 
continues  to  remain  in  phase  A.  The  |emperature  which  was  increasing  for  0  <  t  <  2000  begins  to 
decrease  at  t »  2000  sec;  at  t  3304  sec  the  bar  is  sufficiently  cool  for  phase  M'*’  to  nucleate  and 
begin  to  grow,  until  eventually  at  t  ••  3834  sec  the  entire  bar  consists  of  M**^. 

The  qualitative  features  of!  the  elongation  history  shown  in  Figure  10<a)  can  be 
understood  from  the  preceding  discussion  (Figure  10<b))  by  keeping  in  mind  that  M'  is  the 
low-strain  phase,  A  is  the  intermediate-strain  phase  and  M*^  is  the  high-strain  phase.  During  the 
initial  stages  0  <  t  <  97%  and  978  <  t  <  1026,  when  M'  is  disappearing,  first  due  to  the  growth  of 
M^  and  then  due  to  the  growth  of  both  M*^  and  A,  the  bar  gets  longer.  During  the  stage  1136  <  t 

<  1314  the  bar  is  transforming  front  M"*"  to  A  and  so  the  bar  gets  shorter.  Next,  for  1314  <  t  < 
3304,  the  bar  remain;  in  phase  A  and  so  its  length  does  not  change  appreciably.  Finally,  for  3304 

<  t  <  3834,  the  bar  transforms  from  A  to  M^  and  so  it  gets  longer  again. 

<v).  Finally  we  simulate  two  recent  experiments  carried  out  by  MUller  and  Xu  [26]  on  a 

Cu-Zn-Al  shape  memory  alloy.  Consider  a  bar  of  austenite  at  an  initial  temperature  >  A  .  In  the 

s 

first  simulation,  the  bar  is  subjected  to  an  isothermal  mechanical  loading  during  which  the 
elongation  is  increased  monotonically  until  M'*’  martensite  has  nucleated  and  begun  to  grow; 
then,  before  the  bar  has  transformed  completely  to  martensite,  it  is  unloaded  by  decreasing  the 
elongation  back  to  zero.  The  experiment  is  now  repeated,  with  unloading  commencing  at 
different  instants.  Figure  11(a)  shows  the  result  of  this  simulation  (carried  out  at  6  =  330^  and 

16 1  =  8.333  X 10  m/sec)  and  may  be  compared  with  the  experimental  results  reported  in  Figure 
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13(b)  of  Muller  and  Xu  [26]. 

In  the  second  simulation,  the  bar  is  subjected  to  an  isothermal  mechanical  loading  during 
which  the  elongation  is  increased  until  the  bar  has  transformed  completely  into  M'*’  martensite. 
Next,  the  elongation  is  decreased  monotonically  until  austenite  has  been  nucleated  and  begun  to 
grow;  then,  before  the  bar  has  transformed  completely  back  to  austenite,  it  is  reloaded  by 
increasing  the  elongation.  The  experiment  is  now  repeated,  with  reloading  commencing  at 
different  instances.  Figure  11(b)  shows  the  result  of  this  simulation  which  was  also  carried  out  at 

0  =  330®K  and  I&  =  8.333  x  10"^  m/sec ;  cf.  the  experimental  results  reported  in  Figure  14(b)  of 
MUllerofXu[26]. 
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Appendix:  Restrictions  on  the  material  parameters. 

Here  we  shall  list  all  of  the  inequalities  not  displayed  previously  which  the  material 
parameters  must  satisfy.  According  to  the  statement  below  (3.12)  the  equations  of  the  boundaries 
of  the  regions  Pj  in  the  (y,0)-plane  are  given  by 

^2(0)  =  -  -  Yj.  +  a(0  -  0.J,)  for  0  <  0  <  0j^, 

^2(0)  =  -  Oj^(0)/p  +  a(0  -  aj.)  for  9^  <  B  <  9j^. 

?3(8)  =  Oj^<0)/p  +  a(0  -  02-)  for  <  0  <  0j^. 

^4(0)=  Oj^(0)/jji+  Y'i-+  a(0-0j)  for  O<0<0j^, 

where  the  stress-levels  o^<0)  and  c<^(0)  are  given  by  (3.12).  In  order  that  the  corresponding 
straight  lines  in  the  y.  0*plane  be  ananged  as  shown  in  Figure  1,  it  is  necessary  that  ^4(8)  > 
^3(0)  >  Y2(0)  >Yi(0)>  -1  for  0jj^<  0  <  0j^  and  that  ^4(0)  >  ^^(0)  >  -1  for  O<0<0jjj. 
These  inequalities  can  be  expressed,  upon  usmg  (A.1),  as 

0  <  Oj^(0)  <  Ojj^(0)  +  nY.p  <  p  +  pa(0-0.p)  for 

(A.2) 

0  <  0^(0)  +  py^.  <  p  +  pa(0-Oj)  for  0  <  0  <  0^^. 

Next,  since  we  assumed  in  Section  3  that  all  three  phases  M',  A  and  M'*'  exist  when  o  = 

0  and  0  =  0p  it  is  necessary  that  the  corresponding  strains  y  =  -Yj,  0  and  fj  lie  in  the 

appropriate  strain  ranges  as  defined  by  Figure  1.  In  view  of  (A.1)  and  (A.2),  one  finds  that  this 
holds  if  and  only  if 

yiji  <  1,  ^  f A.3) 

We  turn  finally  to  the  issue  of  extending  the  Helmholtz  free-energy  function  (3.10)  to  the 
unshaded  ("unstable”)  region  of  the  (y,  0)-plane  shown  in  Figure  1.  Even  though  we  do  not  need 
an  expression  for  t|r  on  this  region,  it  is  still  necessary  to  know  that  (3.10)  can  extended  to  that 
region  in  the  manner  previously  assumed  (see  paragraph  below  (3.12)).  The  ability  to  do  this  is 
equivalent  to  the  ability  to  connect  each  adjacent  pair  of  rising  branches  of  the  stress-strain  curve 
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in  Figure  2  by  a  declining  branch  with  prescribed  area  under  it.  Since  the  stress-strain  curve  is  to 
be  declining  for  strains  in  the  intervals  (yj(0),  <y3<0).  Y4<0))  when  0j^<  0  <  0j^,  and 

on  (Yj<0),  Y4(0))  when  0  <  0  <  0^^,  it  is  necessary  that 

Om<0)  >  0^(0)  for  ‘ ^ 

Next,  as  is  readily  seen  from  Figure  2(a).  for  0^<  0<0M.  the  area  under  the  graph  of 
o(«,  0)  between  Y  ?3(B)  and  Y  =  ^4(^)>  Qtust  necessarily  lie  between  the  areas  of  the  two 
rectangles  with  the  same  base  ^^3(0)>  ^4^^)^  with  heights  Oj^(0)  and  (>j^(0)'  A  similar 
restriction  applies  to  the  area  between  Y  =  ?j<0)  and  Y  =  and  for  0  <  0  <  0^^^^  to  the  area 
between  Y  =  ?2(0)  and  y  =  ^4<0)-  Thus  it  is  necessary  that 

o^<0)(V0)-t3(0))  <  PV<Y4(0).0)-PV<V®>*®>  <  Oj^(0)(V®>’V®^)' 
On,<0)(V0>*V®^)  PV<V®>-®>*PV<V^).0)  <  • 

where  the  first  two  sets  of  inequalities  in  <A.S)  hold  for  0qj  <  0  <  0[y{>  while  the  last  set  holds  for 
0  <  0  <  0JJJ.  Conversely,  given  two  points  (^3(0).  Oj^(0))  and  (?4(0),  Oj^(0))  in  the  (y,  o)- 
plane,  with  ^4(0)  >  ^3(0),  a  sujficient  condition  for  the  existence  of  a  continuous  decreasing 
function  o(a,  0)  connecting  these  two  points,  which  is  such  that  the  area  under  it  is  p  v(Y4(0)>  0) 
-  p  V(?3(0).  0)f  is  that  (A.4)^  and  <A.5)2  hold.  The  requirements  (A.4)lr  (A.5)  are  therefore 
necessary  and  sufficient  for  the  extendability  of  the  Helmholtz 
unstable  region. 

The  inequalities  <A.5)  can  be  expressed  equivalently 
Oj^<0)  and  0^(0)  as 


2nYT[o^,(e)  -Oj^(e)] 

"=  2nYT[Oj^(0)-o^(e)] 

•  t*Yx  <  Oni<®>  0 

The  inequalities  (A.2)  >  (A.4)  and  (A.6)  must  be  enforced  on  the  material  model.  They 
can  be  reduced  by  using  (3.12).  (3.13)  into  temperature  independent  inequalities  that  involve 
only  the  material  parameters.  We  shall  not  display  the  resulting  inequalites  here.  These 
inequalities,  as  well  as  (4.3),  are  to  be  imposed  on  the  material  constants  entering  into  our  model. 
One  can  verifiy  that  the  particular  values  (5.6)  of  the  material  constants,  together  with  a  range  of 
values  of  the  four  remaining  parameters  m,  M.  6^  and  «M’  do  satisfy  these  inequalties.  For 

example,  one  possible  set  of  values  of  the  latter  four  parameters  are  m  =  9.7253  x  10*^  /°K, 

M  =  10.1371x10’^  Ac.  =  285°,  0j^=  10,000°;  as  mentioned  previously,  the  particular 
values  of  these  four  material  constants  does  not  affect  the  response  of  the  material. 


for  0^  <  0  <  0^, 
m  M 

for  ®  ^  ®M’ 

for  0  <  0  <  0„. 

m 
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